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Abstract. We show how to compute the probability of any given local con- 
figuration in a random tiling of the plane with dominos. That is, we explicitly 
compute the measures of cylinder sets for the measure of maximal entropy fi 
on the space of tilings of the plane with dominos. 

We construct a measure v on the set of lozenge tilings of the plane, show 
that its entropy is the topological entropy, and compute explicitly the v- 
measures of cylinder sets. 

As applications of these results, we prove that the translation action is 
strongly mixing for fj, and u, and compute the rate of convergence to mixing 
(the correlation between distant events). For the measure v we compute the 
variance of the height function. 

Resume. Soit fi la mesure d'entropie maximale sur l'espace X des pavages du 
plan par des dominos. On calcule explicitement la mesure des sous-ensembles 
cylindriques de X. De meme, on construit une mesure v d'entropie maximale 
sur l'espace X' des pavages du plan par losanges, et on calcule explicitement 
la mesure des sous-ensembles cylindriques. 

Commc application on calcule, pour /i et u, les correlations d'evenements 
distants, ainsi que la v- variance de la fonction "hauteur" sur X 1 . 



1. Introduction 

A domino is a 2 by 1 or a 1 by 2 rectangle, whose vertices have integer coordi- 
nates in the plane. Let X be the space of all tilings of the plane with dominos. Then 
X has a natural topology, where two tilings are close if they agree on a large neigh- 
borhood of the origin. In this topology, X is compact, and Z 2 acts continuously on 
X by translations. 



Burton and Pcmantle |BP] proved that there is a unique invariant measure /x of 
maximal entropy for this action. This measure has the following property: let T 
be any tiling of a finite region R, and Ut the set of tilings in X extending T (we 
call Ut the cylinder set generated by T). Then h(Ut) only depends on R and 
not on the choice T of a tiling of R. 

In this paper we compute ^{Ut) explicitly, for any finite tiling T. We will prove 
in section || that: 

Theorem 1. There is a function P: Z 2 — ► C with the following property. Let E be 
any finite set of disjoint dominos, covering "black" squares b\, . . . ,bk and "white" 
squares W\,... ,Wk- Then (J-(Ue) is the absolute value of the determinant of the 
k x k matrix M = (rriij), where rriij = P(bi — Wj). (Here bi — Wj is the translation 
vector from square Wj to square bi.) 

We call P the coupling function for fi. It can be computed explicitly (see 
Theorem |l6] below), and the values of P are in Q © iQ. 

A similar result holds for tilings of the plane with lozenges. A lozenge is a 
rhombus with side 1, smaller angle tt/3, and vertices in the lattice Z[e 27 ™/ 3 ]. We 

l 
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define a measure v on the space X' of lozenge tilings, which is the limit of uni- 
form measures on periodic tilings, and show that its entropy equals the topological 
entropy of the Z 2 -action. (For a background on measure-theoretic and topological 
entropy see jp|.) We conjecture that v is the unique measure of maximal entropy. 
We compute the coupling function for v in Theorem^: it takes values in Q© 



We give several applications of these two results. First, it is known that the 
translation-action of Z 2 on dominos or lozenges is topologically mixing (if two events 
are far enough apart, one can find a single tiling containing both of them). We show 
here that it is strongly mixing for the measures /i and v, that is, events that are 
far apart are almost independent. In particular we show that the convergence rate 
to independence is at least quadratic in the distance, that is, for any two finite 
tilings Ti and T 2 , and v £ Z 2 , 

fx(U Tl n U (v+T2) ) = v(UtMUt 2 ) + 

and similarly for v; see Theorem [ll]. There is a corresponding lower bound when 
T\ , T2 are single tiles. Correlations for single tiles have been previously computed 



by Fisher and Stephenson [FS| and Stephenson [ptj. 

Another application is to the computation of the variance of the height function 
for v. The height function of a lozenge tiling was introduced by Blote and Hilhorst 
|BHj . We compute for v the variance of the height function in the y-direction 
(Theorem [15]), that is, the variance in height between two vertices separated by 
distance n in the y-direction. This variance is related to the expected number 
of "contour lines" separating the two points (Proposition [l2]). The variance in 
other directions can also be computed with the same methods, although we could 
not obtain a closed form. On the other hand, domino tilings also have a height 
function, but we were not able to compute its variance using the methods here. 

A third application which we will not discuss here is to the computation of 
entropy ent(s,£) of tilings whose height function has a fixed average slope, that 
is, the set of tilings of the plane whose height function ht (-, ■) s atisfies ht(x, y) 



W(0,0) + sx + ty + o(\x\, \y\). This computation is used in CKP ] to give a method 
of counting the approximate number of tilings of any region, by maximizing an 
entropy functional over the region. 

The results of this paper are very general. Using a result of Kasteleyn |K2[ , we 
can compute cylinder set measures (for certain measures of maximal entropy) for 
the space of perfect matchings on any periodic planar graph. Here "periodic" means 
the graph is finite modulo an action of Z 2 by translations. In particular Hurst and 



Green |HG| showed that the Ising model on a planar graph with nearest-neighbor 
interactions can be represented as a matching problem on another planar graph 
(actually they showed this in an equivalent form). Our method therefore allows 
one to compute the local probabilities for the Ising model. Indeed, methods similar 



to those in this paper were used in Montroll, Potts and Ward MPW to compute 
long-range correlations for the Ising model. 

The coupling function which we define is analogous to the Green's function for a 



resistor network (see [BP|). Indeed, subdeterminants of the Green's function matrix 
give probabilities of finding subsets of edges in a random spanning tree. There 
is a well-known connection between domino tilings and spanning trees: Burton 
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and Pemantle [BP| have found, via enumeration of spanning trees, the measures 
of certain cylinder sets for dominos. Our method, which is different, gives the 
measure of all cylinder sets and is adaptable to planar graphs where the spanning 
tree connection is not known to hold. 

The structure of this paper is as follows. Section ^ discusses background: we 
show how to compute the number of tilings as a determinant, and we define /i, v 
and the height function. In section we show how to compute the cylinder set 
measures for finite planar graphs. In section |^ we compute explicitly the coupling 
function for v and give the applications. We chose in this section to concentrate on 
the case of lozenges rather than dominos since lozenges seem to be less prominent in 
the mathematics literature. Section [| restates the results of section ^ for dominos. 

2. Background 

2.1. Matchings on finite planar graphs. Let G be a graph. Two edges of G 
are disjoint if they have no vertices in common. A perfect matching of G is a 
set of disjoint edges which contains all the vertices. If G is finite, connected and 
bipartite, we say G is balanced if there are the same number of vertices in each 
half of the bipartition. 

Let 1? denote the graph whose vertices are 1? and whose edges join every pair of 
vertices of distance 1. It is easy to see that domino tilings of the plane correspond 
bijectively to perfect matchings of 1? . The graph 1? is bipartite: color the vertex 
(x, y) black if x + y = mod 2, otherwise color it white. Note that then a domino 
"covers" exactly one vertex of each color. 

Similarly, lozenge tilings correspond bijectively to perfect matchings of a particu- 
lar bipartite planar graph H, the "honeycomb" grid. The graph H is the 1-skeleton 
of the periodic tiling of the plane with regular hexagons (see Figure [j]) . Rather than 



FIGURE 1. A matching on the graph H and some of the corre- 
sponding lozenges (dotted lines). 
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use lozenges as denned in the introduction, we adopt the following coordinates. Let 
i=§-^andy=§ + H^. The set of vertices of H is V = V U V u where 
Vo = {ax + by | a, b 6 Z} are the "black" vertices, and V\ = 1 + Vo are the "white" 
vertices. There is an edge between every pair of vertices at unit distance. We will 
denote the vertex ax + by by the triple (a, b, 0), and a vertex ax + by + 1 by the 
triple (a, b, 1). Hopefully denoting vertices with ordered triples will not confuse the 
reader. Note also that a lozenge has side V3 in these coordinates. 



Kasteleyn [Kl| and Temperley and Fisher [LF] independently computed the 



number of perfect matchings of an m x n grid (in other words, the number of 



domino tilings of an m x n rectangle). Kasteleyn in [Kl] also computed the number 
of matchings on a toroidal graph in the case of dominos. Later, in [K2], he showed 
how to compute the number of perfect matchings of general planar graphs (thereby 
including the case of lozenges) . His method is to count tilings using Pfafhans. We 
outline here an adaptation of his proof for the case of a simply-connected subgraph 
of H. 

We say that a subgraph of H is simply connected if it is the 1-skeleton of a 
simply connected union of "basic hexagons" of H. The following is a special case 



of the results of [K2 



Theorem 2. Let H' be a finite simply connected subgraph of H . The number of 
perfect matchings of H' equals y/\ det(A)\, where A is the adjacency matrix of H' . 

Here is a sketch of the proof: since the graph H' = (V' , E') is bipartite, we have 
V = V ' U V{ (where V Q ' C V , V{ C V%) and E' C V£ x V{. Order the vertices V so 
that those in Vq comes first. We then can write 



(2.1) 



B 
B l 



where B : Rl y o'l _> Rl^i'l. We can assume that Vq and V{ have the same cardinality 
n (that is, H' is balanced), for otherwise there are no perfect matchings (and det(A) 
is clearly 0). Then the determinant of the square matrix B is 

(2.2) det B = sgn(cr)a lCT(1) • ■ • a na{n) 

a 

where A = (a,j) and the sum is over bijections a: V ' — + V{. Each nonzero term 
in the sum corresponds to a matching. One need only check (and this is the most 
interesting part, although we won't do it here) that each term has the same sign. 
Thus | det(£)| = ^\det(A)\ = # of matchings. 

This same method also works for dominos, except that one must modify the 



adjacency matrix to make the signs in (2.2) work out right. One way to modify it 



is to put weight i = y— 1 on vertical edges [W|. 

It is slightly more complicated to count matchings on a graph embedded on a 
torus. Kasteleyn showed that, for a graph on a torus, one can count matchings 
using a sum of 4 Pfaffians (see Lemma below). 

2.2. The measure of maximal entropy. An important issue not dealt with by 
Kasteleyn is the appropriate sense in which random tilings of large finite regions in 
the plane resemble random tilings of the whole plane. 

Indeed, as the work of [ EKLP| shows, the local configuration of the boundary 



of a region has a drastic effect on the number of perfect matchings, and hence one 
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must take care when computing the measure \x as a limit of measures on tilings of 
finite regions. 

A tiling with free boundary conditions, or simply free tiling, of a region 
R is a tiling of _R, where the tiles are allowed to protrude beyond the boundary of 
R. The topological entropy of domino tilings is by definition 

H top = lim -^loglSfcl, 

k^ca K 

where Sk is the set of free tilings of a k x k square. Here the square regions S k 
can be replaced by any sequence of "sufficiently nice" regions, for example, convex 
regions containing squares of size tending to infinity. In this case one divides by 
the area of the A:-th region, not by k 2 . The same definition also applies to lozenge 
tilings (replace the k x k square with for example the set of vertices in a k x k 
square centered at the origin). 

The entropy H(ijl') of a translation- invariant measure fi' on domino tilings of 
the plane can be defined by 

H( f /) = lim 1 V -»'(U T )lo gf x'(U T ). 

We have H(pf) < H top for any translation-invariant measure // on X, since, as the 
reader may readily show, for a finite set of real numbers {a,i} satisfying a, > and 
53 flt = 1) the quantity — ^3 a i l°g a i is maximized when the are equal. 

For any measure // defined on the finite set X(S) of tilings of a finite region S 1 , 
the entropy (per site) of /i' is 

H (»') = -W\ E log /A*), 

|i| sex(s) 

where f/(s) is the probability of the tiling s. 

2.2.1. The domino measure. Burton and Pemantle showed: 



Theorem 3 ([BP]). There is a unique measure of maximal entropy fi for domino 



tilings of the plane. The entropy of /i equals the topological entropy H top (X). 

Let l? m n be the quotient of I? by the lattice of translations generated by (2m, 0) 
and (0, 2n). Then n is a graph (on a torus) with 4mn vertices. Let y~ m ,n denote 
the uniform measure on perfect matchings of n . 

Let /x* be a weak limit of \i m ,n as to, n — > oo. We claim that fi* = /i. To see 



this, Kasteleyn Kl] showed that the entropy of fi m .n converges as m,n — > oo to 
the entropy of fA. So it suffices to show the limit of the entropies is the entropy of 
the limit. Let Wk by a fc-by-fc window in Z 2 (or on n when to and n are larger 
than k). Let |U TO , n (Wfc) denote the restriction (projection) of ^ m ,n to Wk- Then 
A i m,n(W / fc) converges to n*(Wk) by weak convergence. Since there are only a finite 
number of configurations on Wk , 

lim H^ m , l (W k ))=H( j i*(W k )). 

m,n — >oo 

Since Zf (/i m , n (Wfc)) > H(^ m<n ) we have 

H(fi*(W k ))> lim H(fx m>n ) = Htop- 

m,n — *-oo 

Taking limits over k gives H(ft*) — H t0 p and therefore by the theorem fi* = /i. 



6 



RICHARD KENYON 



We would like to thank Jeff Steif and Jim Propp for the idea behind this argu- 
ment. 

2.2.2. The lozenge measure. Burton and Pemantle's proof of Theorem || above uses 
a connection between spanning trees and domino tilings. Lozenge tilings are in fact 
also in bijection with a set of spanning trees, or rather, directed spanning trees 



( "arborescences" ) on a directed triangular grid, see [PW|. However it is not evident 



how to adapt the proof of |BP] to this case. We will show nevertheless that the 
uniform measures on certain toroidal graphs converge to a measure is, and their 
entropy converges to the topological entropy. By the argument following Theorem 
||, the entropy of v equals the topological entropy. The only fact missing is the 
uniqueness of v. 

Let H m n be the quotient of H by the lattice of translations generated by mx 
and ny. Then H m ^ n is a graph with 2mn vertices. Let is m ,n be the uniform measure 
on matchings of H m ^ n . 

Lemma 4. The entropy of v n ^ n converges to the topological entropy of lozenge 
tilings of the plane. 

Proof. The proof is a simple variant of an unpublished trick due to G. Kuperberg. 
Let F n be the set of free lozenge tilings of the region R n , where R n consists of an 
equilateral triangle of side n^/i (recall that a lozenge has side length . Then F n 
is partitioned into equivalence classes, two tilings being equivalent if they have the 
same set of tiles which protrude beyond the boundary. The number of equivalence 
classes is less than C n for some constant C, so some equivalence class B n must 
satisfy 

C n \B n \ > \F n \. 

Reflect the region R n along one of its edges c. Given any two tilings Ti, T% e B ni 
there is a tiling of the union of R n and its reflection, where the R n is filled with T\ 
and its reflection is filled with the reflection of T^. This is because the tiles which 
protrude across the boundary edge c are fixed under this reflection. 

Take the union of 18 reflections of R n so as to make the parallelogram of Figure |[ 
Given any 18 elements of B n , we can fill them in the 18 triangles to make a tiling 
of this whole parallelogram. Furthermore, the configurations of tiles protruding 
beyond opposite boundaries of this parallelogram are translates of each other. As 
a consequence we can glue opposite sides of the parallelogram together to make a 
tiling of the 3n by 3n torus (tilings of which correspond to perfect matchings of 
H^n^n)- We therefore have (where N 3n is the set of tilings of H 3n 3n ) 

(C- n \F n \) 18 < \B n \ 18 < \N 3n \ < \F n \ 18 . 

This implies that 

lim \ log \F n \= lim -i^loglAq, 

which completes the proof. □ 

We will see later that H{v m ^ n ) and H{v n ^ n ) have the same limit as m,n — > oo, 
so the restriction to the tori H 3rit3n is unimportant. In fact we will show in section 
4.4 that there is a unique weak limit to the v m ,n- Let v be this limit. We make the 



following 

Conjecture. The measure v is the unique measure of maximal entropy for lozenge 
tilings of the plane. 
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Figure 2. The 18 triangles make a torus. 



2.3. Height functions. For any perfect matching of a simply connected subgraph 
H' of H one can associate an integer- valued height function to the faces of H' . 
(By face of H' we mean a basic hexagon which contains an edge of H' .) 

Let T be a perfect matching of H' . Pick a face xq arbitrarily and assign it height 
h(xo) = 0. For any other face x, take a path xq,xi, . . . , x n = x of faces, where for 
each j, Xj + \ is adjacent to Xj along an unmatched edge. Let w = e 27 ™/ 3 . Then 
h(xj + i) = h{xj) ± 1, where h increases by 1 if Xj+\ — Xj is in directions i,iw, or 
iw 2 and h decreases by 1 if Xj+i — Xj is in directions — i, —iw or —iw 2 . This defines 
a height on each face, which can easily be shown to be independent of the choice of 
path to that face. Figure |3| shows a height function for the matching in Figure [l]. 

For dominos there is a similar definition; see Q . 

For a simply connected subgraph H ' which has a perfect matching, the height 
function along the faces bounding the region (that is, faces outside H' but contain- 
ing an edge of H') is well defined up to an additive constant and is independent 
of the matching Q (rather, this is true if the union of H and these faces is still 
simply connected). 



3. The measure of a cylinder set 

The results in this section are stated for lozenges, but they apply verbatim for 
dominos, on condition that one puts weight i = \/— T on vertical edges in the 
adjacency matrix B. 

Let H' be a finite, balanced, simply connected subgraph of H, with vertices 
{vi, . . . ,v n } € Vq and {v[, . . . ,v' n } G V\. Let E = {ei, . . . ,ek} be a subset of 
disjoint edges of H', with ej = v Pj v' q . for j = 1, . . . , k. 

In the formula ( |2.2|) , each nonzero term in the sum corresponds to a perfect 
matching; those terms whose matchings contain all the edges in E are precisely the 
terms which contain the subproduct 



def 

a E — a pigi a P2q2 ' ' ' a Pkqk ■ 
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Figure 3. Height function for lozenges. 

The sum of the terms containing this product is (up to sign) the product of ag 
with the determinant of the cofactor Be, that is, the determinant of the matrix 
obtained from B by removing all rows pi, ■ ■ ■ ,Pk and all columns qi,... ,qu- This 
can be proved inductively as follows. Expand det(B) along the row p^: 

det(B) = (-l) Pl+1 (a Pl i detB pul - a Pl2 det B P1 . 2 + ... ± a pin det B Pu „). 

The only term that contains a piqi is the term for column q±, which is the term 
(— l) Pl+9l a Pl9l det-B Pl9l . Now expand B piqi along row p 2 , and so on. Thus the 
sum of the terms containing oe is (— 1)^ Pj+qj a,E det(_Be). 
Since \oe\ = 1 this proves 

Proposition 5. The number of perfect matchings containing all edges in E is 

|det(S B )|. 

Later on the sign of det(B^) will be important to us. We note therefore that 
the number of perfect matchings containing all the edges in E is 

(_X)£(w+«*) aB det(S B )sign(detB). 

In a more general setting (e.g. in flCKPfl ) one might wish to consider weighting 
the edges non-trivially, in which case the formula \cle det(S^)| counts each matching 
according to the product of the weights on its edges. 

Note that the quantity det(i?E) does not depend on the actual set of edges in E, 
but only on the set of vertices involved. This shows, as expected, that each matching 
of the vertices of E can be extended to the same number of perfect matchings. 

The matrix Be is an (n — k) x (n — k) cofactor of B; its determinant is equal to 
det(_B) times the determinant of the k x k cofactor of the inverse of B, (B~ 1 )e*, 
where E* is the set of rows and columns not involved in E. So we have 
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Theorem 6. The number of perfect matchings containing all edges in E is 

|det((B- 1 ) £ »)det(B)|. 

That is, 

Hh>(Ue) = |det((S- 1 ) E ,)l- 

More precisely, 

»w(U E ) = (-l)E« + «a B det((B-V)- 



The advantage of this over Proposition ^| is that the computation is a determinant 
of size k rather than n — k. In particular if E consists of a single edge ViVj, the 
probability that that edge is in a random matching is the absolute value of the ij-ih 
entry of B~ x , 

The rest of the paper consists of the computation of B~ l and applications of 
Theorem |[ 

4. Lozenge tilings 

4.1. The torus. Recall that H m ^ n is the toroidal graph H/(mZx + nLy), and the 
vertex rx + sy + t is designated by the triple (r, s, t) S Z/mZ x Z/nZ x {0, 1}. Let 
A\ be the adjacency matrix of H m n . Let A-i be the matrix obtained from A\ by 
changing the sign of the entries corresponding to edges from (m— 1, k, 1) to (0, k, 0) 
for each k £ Z/nZ. Let ^3 be the matrix obtained from A\ by changing the sign of 
the entries corresponding to edges from (k, n— 1, 1) to (fc, 0, 0) for each k S Z/mZ. 
Let A4 be the matrix obtained from A\ by changing the sign of both these sets of 
entries. 

Let Bi, B2, B3, B4 be the matrices obtained from A%, A2, A3, A4, respectively, as 



in (2.1) 



Following Kasteleyn, on condition that m, n are both even, we have 
Lemma 7 (|K1|). The number of perfect matchings Z m . n of H m ,n is given by 



(4.1) Z m , n = -(- dct B 1 + dct B 2 + det B 3 + dct B 4 ). 

A term for each perfect matching appears in each of these four determinants. 
The signs are arranged so that each perfect matching is counted three times with 



sign + and once with sign — (see Kl]). For other parities of m and n a similar 
equation holds but with different signs. Throughout the rest of this section we will 
assume m and n are both even and neither is divisible by 3. This is simply a matter 
of convenience; the computations go through in the other cases but require slightly 
different arguments. 

Theorem ^| also applies to tilings of the torus. The quantity 

detBi det^f 1 )^ 

contains a term for each matching containing all edges of E. The sign of these 
terms is the same as the sign of the corresponding terms in det(Bi) (this easily 
follows from the proof). This is also true for #2,-83,-84, and so 
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Lemma 8. The number of perfect matchings of T rnn containing a given set of 
disjoint edges E is given by 

(4.2) ^ l-det-BidetOBf 1 )*;. + det B 2 det(-B^%.+ 

+ dctB 3 dct(B 3 " 1 ) E . +dct J B 4 dct(B 4 : 1 ) i j»| . 



4.2. Computation of the eigensystems. The eigenvectors of the Aj are func- 
tions 

f(r, s,t) — c t z r w s , and /'(r, s, t) = c t z r w s 

where 

c = vT+ + 

Cl = Vl + 2 + w 
c'o = c 
ci = -ci 

(for either choice of sign of the square root), and 





= 1, 


W n 


= 1 


for Ai 


z' m 


= -1, 


w n 


= 1 


for A 2 


z m 


= 1, 


w n 


= -1 


for A 3 


z ra 


= -1, 


w n 


= -1 


for A 4 



The corresponding eigenvalues are 

A = c ci 
A' = — c ci. 

Consider for example the matrix A\. Let 







W\ 




Ck.tfi 




Ck,e,i 


= y/l + z^ + wi 


c 'k,e,o 




c 'k,t,i 


= -Ck,e,i 



(any choice of sign for the square roots will do). Then fk,e(r,s,t) = Ck,e,tZi r Wi s 
and f' ke {r, s, t) = c' ke t Zi r w[ s are eigenfunctions. The functions fk,e + f' k t are zero 
on Vi, and run through a basis for all functions on Vq (the exponentials) as (k,£) 
runs through {0, . . . , m — 1} x {0, . . . , n — 1}. Similarly the functions fk.i — f k £ 
are zero on Vq and run through a basis for functions on V\ . Thus the collection of 
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fk,e an d f'k i is a basis of eigenvectors for the matrices A\. In particular we have 

m— 1 n— 1 

det^O = 1] II a M*m 

k=0 1=0 
ra— 1 n— 1 

k=0 £=0 
m— 1 n— 1 

(4.3) = ]J II ( X + e 27rtk/rn + e 27 ™ £ /") 2 

k=0 £=0 

(we removed the initial minus sign in this last equality since there are an even 
number of terms, and also in the second product we replaced k,£ with — k, —t). 
Similarly one can compute 

det(A 2 ) = [J(l + e ^(2fc+i)/™ + e 2mt / n ){l + e -^k+i)/ m + e -2*u/ n ^ 
k,e 

(4.4) = + e 7ri ( 2k + 1 )/ m + e 27ri£/ri-j2 

k,e 

and furthermore 

(4.5) det(A 3 ) = JJ(1 + e 2»ifc/m + e? ri(2/+l)/nj2 

(4.6) det(Ai) = ]J(i + e ^( 2fe + 1 )/™ + e ^( 2£ + 1 )/«) 2 . 

These four determinants are non-zero by our assumption that m and n are nonzero 
modulo 3 (three complex numbers of modulus 1 sum to zero only if they are a 
rotation of the set {1, e 27 ™/ 3 , e 47 ™/ 3 }). 

4.3. Computation of the inverses of B\, . . . , B4. It suffices to invert the A,. 
Define the vector 

1 if (x,y,t) = (0,0,0) 
otherwise. 



5o,o,o{x,y,t) 
aatrix A\. W 



Consider first the matrix Ai. We can write 

fk,e{x,y,t) + f k Ax,y,t) 



2mn ^ c fc ,£,o 



Thus the vector 



(4.7) P (a; )y ,t) = J_V— [ 



f k ,e{x,y,t) fk tt (x,y,t) 



A/t.£ A', 



has the property that 

^i-Po,o,o(a;,2/,*) = ^0,0,0(2;, J/,*)- 
The preimages of delta functions at other vertices can be similarly defined: define 
Pr,sfi{x,y,t) = Poflfi{x - r,y - s,t); then 

1 if (x, y, t) = (r, s, 0) 



^i-Pr,s,o(a;,y,*) = S r>3t0 (x,y,t) 



otherwise. 



Thus the function P r , s ,o gives the entries of the (r, s, 0)— column of the matrix A 1 1 . 
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For notational simplicity we will write P^'(x,y,t) = Pofi t o(x,y,t) (the super- 
script (1) refers to the fact that it comes from matrix A\\ note that pW also 
depends on m and n). 

If we simplify ( ff.Tj ) above: 

1 ( fk,t{x,y,t) - f ke (x,y,t)' 



pW(x,y,t) 



2mn ^— ' c k ,e,o 

k.l ' ' 



A 



kj 



we see that P^(x,y, 0) = for all x, y. On the other hand plugging in for 

/, c, A, Zi,wi we have 



(4.8) 



pW( X ,y,l) = J-J2 T 

mn L — ' 1 

k,t 



j j D — 2irik/rn _|_ ^ — 2'Kiijn ' 



Similar expressions hold for the inverses of A2, A3, A^, except that the sum is 
over a different set of angles. For example for Ai we have 



mn L — ' 1 



^Ki{2k+l)x / 'm ^2-Kiiy / 'n 



jjljl t ^ _|_ g— 7ri(2fc+l)/m _)_ p—2-Kit/n ' 



4.4. The limiting measure. Here we compute the limit of the ratio of (4.2) and 
(hi), that is, the limit of the cylinder set measures v m>n {UT)- We still assume 
m, n ^ mod 3 so that the determinants in ( [4.3| )-( |i~6| ) are nonzero. 

We will show below that the quantities det((B^ 1 )£;* ) for j = 1,2,3,4 all converge 
to the same value as m,n — > 00. Then (if we remove the absolute value sign in 



(1-2)) the ratio of (4.2) and (4.1) is a weighted average of these four quantities, 
with weights ± det Bj/Z m ^ n . These weights are all in the interval [—1,1] since 
Z m ,n > |detBj| for each j: recall that Z m>n and |det-B,-| count the same objects 
except that some signs in |det-Bj| are negative. Since the weights sum to 1, the 
weighted average is also converges to the same value as each det((P^ 1 )£;«). 

It suffices to show that, for fixed x, y, for each j — 1,2,3,4 the quantities 
PW(x, y, 1) tend to the same value (recall that the matrix B~ l has entries which 
are the values of P")). 

For a fixed integers x, y, the function 



-i<f> 



is continuous on ([0,2tt] x [0,2tt]) \ {(2tt/3, 4tt/3) U (4tt/3, 2tt/3)}, and has a pole 
at the two points (2n/3, 47r/3) and (in/3, 2n/3). The expression (l.g) and the 
corresponding expressions for the other P^ are each Riemann sums for the integral 
of p over [0, 2tt] x [0, 2ir}. Fix 5 > and let Ns be a S x (5-neighborhood of the pole 
(#Ch<A)) = (27r/3, 47r/3). Outside Ns (and the corresponding neighborhood of the 
other pole) the Riemann sums converge to the appropriate integral. We must show 
that on Ns the sums are small (in the sense that they tends to zero with S as 
m, n — > 00). 

On Ns, we use the Taylor expansion of the denominator: 



1 

mn 



(e,<p)eN 5 



< 
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— mn ^N s |- 4e -'fo(0-0 o )—ie- i *o(0_0 o )+O((6l-eo) 2 ,(0-0o) 2 )| 



JL V 



1 



1 

mn 



mn ^N s | e -2W3(0_0 o ) +e -4W3(0_0 o ) +O ((e_e o )2 ;( 0_0 o )2)| 

T,N S i e -2^/ 3{e ^ eo) l e -4^ /3(<f> _ Ml \ + 0(0 - 00, ^ - <M- 
(The reason we can take the big-O term out of the denominator is that 

| e -2W3(0 _ + e -4«/3(^ _ 0O )| > co71si (| _ Q | + |^ _ Q |) 

for 1 — 6q I and | <^> — 4>q | small; we can then move it out of the summation altogether 
because of the factor -^-.) 

mn ' 

Since m, n ^ mod 3, we have, in the sum for each P^) , that 9 — 9q is of the 
form 9 — 9q = 7T (-^ — |) = Kj— and similarly (f> — 4>o = where A;' and £' are not 
multiples of 3. Thus the sums are each bounded by one of the form 



77.77. 51 |_fc _l X 4™/3| + 



, : ,-e 4 W 3 | 

I om on l 



where the sum is over 



(M)ez 2 | (^-,^-)€JV 5 \{(0,0)} 

om on 

This sum is easily shown to be bounded independently of m, n and tending to 
zero as S — > 0: this follows from the integrability near the origin of the function 
\ x+e <Li/3 y \ w hh respect to dxdy. 

Theorem 9. For a finite matching E covering black vertices {bi, . . . , bk} and white 
vertices {w\, .. . ,Wk\, we have v{Ue) = | detP(iOj — bj)\, where 



i r n r n ( 



-id _|_ g— itf> ' 

Proof. This follows from Theorem ^ and the convergence argument above. □ 



4.5. The coupling function for lozenge tilings of the plane. We call P(x, y, t) 
the coupling function for lozenge tilings of the plane. Let z — e~~ %e and w — e -4 *. 
Then we obtain the contour integral 



(4.9) 



1 

47^ 



dz dw 



s i xS i (1 + z + w) (-iz) (-iw) 
The function P has all the symmetries of the graph H , that is, 



(4.10) 



P(x,y,l) 



P(- 
P(- 



-y 
y - 



1,2,1) 

i,y,i) = 



■-P(y,- 

P(x,- 



-x ■ 
■x - 



■y 
y- 



i,i). 



These are obtained by noticing that (4.9) is invariant under interchanging z and u>, 



, zw 



1 ), or under combinations of these. 



or under the substitution (z,w) — > («T 

We can explicitly evaluate the integral (4.9) as follows. This is slightly easier to 
calculate in the case x < —1. The other values can then be obtained by (4.1C). If 
we fix w and integrate over z, there is a pole inside the unit circle when |1 + w\ < 1, 
that is, when — e (2tt/3, 47r/3). The residue of z^ x ^ 1 /(l + z + w) at this pole is 
(— 1— w)^ x ^ 1 . For —(f) ^ (27r/3, 47t/3) there is no pole inside the circle and therefore 
the integral is zero. 
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So when x < — 1, 



Pfo j/, !) = —(-!)-*- 1 



(4.11) 

For example when x = — 1 we have 



(1 + wy^dw. 



(4.12) 



^(-l >tf ,l) = _ 



,2TTiy/Z _ A-niy/3 



,V3 



where Cj, = 0, 1, — 1 if y = 0, — 1, 1 mod 3. Some values of P(x,y, 1) are plotted 
in Figure [|. To compute these values, start from ( 4.12Q , which, along with the 




.1-1/3 



^3 72+1/3 



, 3T/2-1/3 



-2T+2/3 



Figure 4. Values of P(x,y,t). Here r = y/3/(2n) 



symmetries ( 4. 10|) gives the value on 3 lines, and use the fact that the sum of the 
three values adjacent to any vertex (except the origin) is zero (since AP = at 
these vertices). For example the equation 

P(-fc, 0, 1) + P(-k + 1,-1,1) + P(-k + 1, 0, 1) = 

determines P{—k, 0, 1) for k > inductively for increasing k (note as before that 
P(-fe,-l,l) = P(-l ) -fc,l) = ^). 

4.6. Some local probabilities for lozenges. From Theorem ^ and equation 
(1.11) with x = —1,J/ = 0, the probability of the edge (0, 0, 0)(0, 0, 1) being in 
a random matching is given by P(0, 0,1) = ^, as expected by symmetry. The 
probability of both the edges (0, 0, 0)(0, 0, 1) and (0, 1, 0)(0, 1, 1) occurring is given 

by 

P ,o.o(0,0,l) P ,o.o(0,l,l) 
Po,i >o (0,0,l) P Oll ,o(0,l,l) 
The probability of the three edges 

(0,0,0) (0,0,1), (0,1,0)(-1,1,1) and (-1, 1, 0)(-l, 0, 1) 



x/3 

2-7T , 



9 07T 



LOCAL STATISTICS OF LATTICE DIMERS 



15 



(which bound a hexagon) all occurring, is: 





1 


1 


2n 


3 


3 


1 


Va 


1 


3 


2?r 


3 


1 


1 


V3 


3 


3 


2tt 



2 

27 



V3 
Gtt 



3V3 
8^3" 



14.5%. 



(Such a hexag on co rresponds to a local minimum for the height function.) 

Note from (|4.12[) that _P(-3fc,3fc,l) = P(-l,-3fc,l) = 0. As a consequence, 
the set of edges E = {e n } ne z where e„ is the edge (— 3n, 3n, 0)(— 3n, 3n, 1) is an 
"independent" set: the probability of any fixed subset of size k occurring in a 
random matching is exactly (l/3) fc . Is there any simple explanation for this fact? 

We can also compute long-range correlations. The probability of the horizontal 
edge (0, 0, 0)(0, 0, 1) and the horizontal edge (— n, n, 0)(— n, n, 1) both occurring is 
given by the determinant 

1 



P(n,-n,l) 



P(-n,n,l) 
i 



9 



P{-n,n,iy 



By Q4.10| ) and ( [4.12[ ), P(-n, n, 1) = P(-l, -n, 1) = c„V3/(27rn). So the probability 
of these two edges both occurring tends quadratically (but not faster) to 
generally we have 



More 



Proposition 10. 



O 



1 



P(x,y,l 

Proof. We show for x < — 1 that P(x, y, t) 
try of P the result will then follow. From flt.ll ) we have 

|P(^y,i)l 



0(|iy); by the dihedral 6- fold symme- 



< 



1 

27 
2 

2^ 
1 



< - 



i/3 



-\-w\ x ~ 1 \dw\ 
(2cos(0/2))- x - 1 d9 



2tt/3 



r 



y/1 WW 
1 



-.dt 



where we used the substitution t — 2 008(0/2). 



□ 



(4.13) 



Let Ei and E^ be two finite collections of edges. The probability that edges 
E\ U E2 are present in a matching is given by a determinant 

Ci D 1 
D 2 C 2 

where det(Ci) = ^(Ue^ and det(C2) = [i(Ue 2 )- The entries of D\ and D 2 are the 
values of P connecting points of E\ to E 2 - If the distance between E\ and E 2 is 
> n then these entries are all O(^). Let J] e(<r)a 1(T ( 1 ) . . . a kcr ^ be the expansion 
of the determinant ( 4.13| ). Suppose a ia ^ is an entry from the submatrix D±, that 
is, i < \Ei\ and a(i) > \E\\. Then there must exist an index j with j > \E\\ 
and a(j) < \Ei\ (since a is a bijection), that is, ijcrlj) is an entry of D 2 . We may 
similarly conclude that each term in the sum for the determinant has the same 
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number of entries from D\ as from £>2- Therefore the determinant (4.13) is equal 
todet(Ci)det(C7 2 ) + 0(^). 
This proves 

Theorem 11. If the distance between E± and Ei is at least n, then 

n{U El n Ue 2 ) = fJ.(U El )n(U E2 ) + O(-o)- 

Note that the comments before Proposition ([l0|) give a corresponding quadratic 
lower bound if E\ and Ei are single horizontal tiles separated by — (3n + l)x + 
(3n + l)y. 

An argument similar to Theorem pT| shows that v is mixing of all orders, that 
is, given k events E±, . . . ,Ek each of which is separated from the others, the joint 
probability converges to the product of the individual probabilities as the distance 
between them tends to infinity. 

4.7. The height distribution. Let h n be the random variable which gives the 
height difference between the face at the origin (by which we mean, the face above 
and right of the point (0,0,0)) and the face at the point (— n, n, 0) in a random 
lozenge tiling of the plane. In this section we compute the variance of h n . 

Since the height is only defined up to an additive constant, we may as well 
assume that the height on the face at the origin is zero. 

One motivation for computing this variance is that it is related to the number 
of cycles in the union of two random matchings as follows. Given two random 
matchings Mi and Ma, their union is a set of disjoint cycles. The difference of their 
height functions hS 1 ' — has the property that it is constant on the connected 
components of the complement of these cycles, and changes by 3 when crossing 
any of these cycles. The height difference across a cycle increases or decreases 
by 3 depending on which of the two possible matchings of the cycle occurs in 
Mi (with the complementary matching necessarily occurring in M2). Since the 
two matchings were chosen randomly, across each cycle the difference of heights 
increases or decreases by 3 independently of what happens at the other cycles. 
(Among all pairs of matchings whose union has this particular cycle structure, each 
matching is equally likely to have either "half of each cycle.) 

Thus for a face / separated from the face at the origin by k cycles, the variance 
in the height difference (that is, the variance a 2 {Ah{f)) = a 2 (h^ (/) - 7i (2) (/))) is 
the same as the variance of a sum of k independent fair coin flips of values ±3. This 
variance is 9k. Since M\ and M2 were chosen randomly, the number y of cycles 
separating / from the origin is also a random variable. We have a 2 (Ah\y) = 9y, 
and 

a 2 (A%) = E((A/ l ) 2 |y) - E(A%) 2 = E((A^) 2 |2/) 
since E(Ah\y) — 0. Furthermore, 

a 2 (Ah) = E((A/i) 2 ) = ^P( y )E((A^) 2 | y ) = ^P(y)9j, = 9E(y). 
y y 
Thus in general the variance of the height differences at a fixed face (for two match- 
ings) is 9 times the expected number of cycles separating that face from the origin. 

Since the two matchings Mi , M2 are chosen independently, the variance of their 
height difference at / is exactly twice the variance of the height of a single matching 
at / (note that the height at / has mean value zero, which follows from the fact 
that the probability of any edge being present is 1/3). We conclude 
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Proposition 12. In the union of two random matchings, the expected number of 
cycles separating two faces /i,/2 is 2s/9, where s is the variance in the height 
(between fi and fi) of a single random matching. 

Let E n denote the set of horizontal edges (— k, k, 0)(— k, k, 1) for k = 1, . . . , n. 

Lemma 13. The value of h n is equal to n — 3r ra , where r n is the number of hori- 
zontal edges from E n in the matching. 

Proof. This follows from the definition of the height function. On the vertical path 
of faces from (0,0,0) to (— n, n, 0), the height increases by 1 for edge of E n not 
present in the matching, and decreases by 2 for every edge present in the matching. 
So the height difference is (n — r n ) — 2r n . □ 



Let M = M n be the n x n matrix M — (rrikj) where 

m kd =P(-\k-j\ ) \k-j\,l). 

Then as we saw, the determinant of M is the probability that all the edges from 
E n are present in a random matching. By symmetry (4.10), when I > 0, 

ceV% 



P{-i,£,l) 



P(-l,-i,l) 
_ V5\ 



2ir£ 



For example when n = 5 we have (with r = 2rr 



Mr 



/ 


1 
3 


— T 


T 

2 





— T 

4 


\ 




— r 


1 

3 


— T 


r 

2 









T 

2 


— T 


1 

3 


— T 


T 

2 









T 

2 


— T 


1 

3 


— T 




V 


— r 
4 





r 
2 


— T 


1 

3 


/ 



For each subset S C E n of cardinality fc, let Ms be the submatrix of M whose 
rows and columns are those of M indexed by S. Then (-1)^ Pl+q > det(Af s ) is the 
probability that the edges corresponding to S are present in a random tiling, where 
the sum is over S. Since pj — qj in our ordering, the probability in question is 
exactly det (Ms)- 

Let 

p( z ) = z n - axz 71 - 1 + a 2 z n - 2 -... + (-l) n a n 

be the characteristic polynomial of M. Then otk is the sum over all subsets S C E n 
of size k of det Mg, that is, the sum over S of the probability that S is in the 
random matching. 

Let q(z) = J2k=o @kZ k where (3k is the probability that there are exactly k edges 
from E n present. A formula of Ch. Jordan [Q relates the au to the f3k'- 



Pk = Q'fc 



jfe + 1 
k 



ctk+i 



k + 2 
k 



Oik+2 



From this we can easily derive the following: 
Theorem 14. We have 

q{z) = (1 - z)>( 



1 



□ 
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Now if p(z) = n"=i( z _ tnen 

n 

q(z)=l[(l-K + X l z). 



»=l 

The expected value of r„ is ^ ^/3fc = ?'(!)■ Wc find 

" A 

so 

q'{\) = q(l) K = Xi = trace ( M ) = ™/ 3 - 
This is as expected since the probability of each edge (— k, k, 0)(— k, k, 1) being 
present is 1/3. 

Similarly, the variance of r n is given by 

a 2 (r„) = £ fc 2 /? fc - (£ k ^ = W> + - ^ 

We have 

= 9(Z) § (l-A J + A 4 zj(l-A,+A,z)' 

and so 

g"(l) = ^AiAj = trace(M) 2 - trace(M 2 ) = (n/3) 2 - tracc(M 2 ). 
To compute trace(M 2 ), note that M is symmetric and the ith row for 1 < i < n 

is: 

/ CjT T _ 1 _ T C»_jT 

^-1'""'2' T '3' T ' 2'"" ' n-i >' 
The inner product of this row with itself is 

I + t 2 (1 + 4 + 4 + 4 + --- + , C "~\J 
9 2 2 4 2 5 2 (n-z) 2 

2/, 1 1 1 C 2 , 

Similar sums (with ^ replaced with 1/3) hold for the rows i = 1 and i = n. 
Summing these up for i = 1, . . . , n yields 

trace(M 2 ) = 2 + r 2 ( 2(n - 1) + 2^ + 2^ + . . . + 2^ 
(4.14) 

= S+AfMi+4+^+-"+7^)-2(i+5+z+-"+ h! 4))- 

9 47r 2 \ 2 2 4 2 (n-1) 2 2 4 n-1 1 
Note that 

1 1 IcJ 

11 1 \ 1 / 1 1 



2 2 3 2 r 2 y 9 V 2 2 3 2 [r/3\< 
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and similarly 



1 1 

1+ 2 + 4 



IfVl 
r 



[r/3\ 



log(r) + 0(l)--(log(r/3)+0(l)) 
Hlog(r) + O(l). 



Plugging this into ( 4.14 ) gives 

trace(Af 2 ) 



n 
3 



log(rc) 



0(1), 



which yields 



Finally we find 



<r 2 (r n ) 



log(n) 



0(1). 



Theorem 15. Let h n be the height difference between the origin and the point 
(— n, n,0). The expected value of h n is 0. The variance in h n is given by a 2 (h n ) = 

+o(i). 



The only previously known result in this vein is in [CEP , where Cohn, Elkies, 
and Propp found a 0(\/n) upper bound for the height variance on a bounded region. 

5. Computation of the coupling function for dominos 

Proof of Theorem ^|. The proof follows from Theorem |6|, the fact that /U m ,„ — > ji 
(see the comments after Theorem [}]), and a convergence argument analogous to 
that preceding Theorem 0. □ 

In fact a computation analogous to the case of lozenges yields: 
Theorem 16. The coupling function for dominos is 



P{x,y) = 



1 



2-7T r 2ir 



J{x8+y<p) 



4t 2 Jo Jo 2 cos(0) + 2i cos(</>) ' 



This statement depends on our choice of weighting the vertical edges of 1? with 
weight i. 

As for lozenges, P(x,y) is zero for black vertices (i.e. when x + y = mod 2). 
Some values are shown in Figure ^| (where the origin is in the lower left corner). 
One can explicitly evaluate, for x > 1, 



P(2x + 1,2a;) 



Tf 



1 1 



1 1 



+ (-1) 



x+l 



2x-l 



and 



P(2x, 2x - 1) = -iP(2as + 1, 2x). 

From these values (and the diagonal symmetry) the others can be obtained induc- 
tively using AP = except at the origin. 
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4 7C 



f _ 5 

71 4 



v 2 3^ 



■i 1 K 



11 
4 71 



11 
71 4 







v -5 4 . 
" 1( 4 + 5b } 



(0,0) 



2_3 
7C 4 



13 
4 



Figure 5. Values of P(x,y) for dominos. 



As an example, the probability of the two edges (0, 0)(1, 0) and (0, 1)(1, 1) both 
occurring is 

I =i 1 

4. 4 = _ 

— j i o 

— 4 8 

The probability of the two edges (0, 0)(1,0) and (0, 1)(0, 2) both occurring is 
given by 

1 



47T 



This probability was also computed in [BP], although their stated value is incorrect 
(they have acknowledged the mistake in their arithmetic). 

Proposition [ic| and Theorem [ll] also apply to dominos. However the computation 
of the variance in the height function cannot be done within this framework since 



there is no simple analogue of Lemma 13 
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